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ABSTRACT:  Oscillation  of  fuel  composition  in  ducted  air  flows  with  oscillatory  injection  of  gaseous 
fuel  was  studied  using  tempentuie  difference  as  scalar  in  flows  involving  radial  injeebon  of  an  oscillatory 
fuel  jet  and  radial  and  axial  injection  of  a  steady  air-fud  jet  with  oscillatory  fuel  composition  into  a  ducted 
air  stream.  Local  measurements  of  mean  and  osdllatory  velocity  and  scalv  in  turbulent  flow  were  carried 
out  with  flow  rates,  frequency  and  amplitude  of  imposed  oscillations  as  variables.  The  decay  of  the 
oscillation  of  the  scalar  was  faster  than  that  ot  mean  velocity  and  mean  scalar,  and  attenuauon  was 
strongest  in  the  shear  layer.  Flow  past  a  bend  attenuated  oscillabons  by  around  40  %  and  past  a  bluff- 
body  downstream  of  a  jet  resulted  in  faster  decay  of  oscillations  than  in  flows  without  a  bluff-body.  The 
influences  of  frequency  in  the  range  between  100  and  220  Hz  and  mass  flow  ratio  of  cross  flow  to  jet 
between  5  and  15  on  attenuabon  of  oscUlatiocts  was  small. 

1.  iKtroducdoa 

Periodic  variaii''n  in  the  composibon  of  an  air-fuel  mixture,  although  not  always  desirable  in  eombusbon 
equipment,  has  useful  applications  as,  for  example,  in  acbve  control  of  combusbon-driven  oscillabons  in 
premixed  flames  [1-4],  O^tUabon  of  the  fuel  composition  so  Uiat  Ute  oscillatory  heat  release  is  out  of  phase 
with  die  combusbon-driven  oscillabons  results  in  their  attenuabon  and  may  be  achieved  by  osdUatory  injection 
of  part  of  the  fuel  into  the  combustor  duct  or  into  a  secondary  air  stream  entering  the  combustor  duct  upstream 
of  the  flame-holder  (3,4],  Effeebve  control  dqiends  on  the  amplitude  and  coherence  of  the  oscillatory  heat  input 
and  the  present  study  is  concerned  with  the  attenuation  of  oscillation  of  fuel  in  flows  of  relevance  to  acbve 
control.  The  oscillabon  of  fuel  composibon  is  quanttfled  in  flows  comprising  an  oscillatory  jet  of  gaseous  fuel 
injected  radially  into  a  Utbe  carrying  secondary  air  and  radial  and  axial  injeebon  of  a  nearly  steady  secondary  air- 
fuel  mixture  with  oscillabng  fuel  composibon  into  the  main  duct  carrying  air.  The  effects  of  a  bend  in  the 
secondary  air  tube  and  a  bluff-body  located  downstream  of  the  secondary  air-fuel  jet  are  also  quanufled.  Fuel 
concentrabon  was  replaced  by  temperature  difference  as  scalar  in  the  measurements  which  include  the  mean 
velocity  and  scalar  and  their  fluctuabon  at  the  imposed  frequency.  The  absolute  tempentiae  difference  between 
the  fluid  streams  was  around  10  %  so  that  the  effect  of  variabon  in  density  was  negligible. 

2 .  The  Eiperimcntat  Set-up  and  Measurement  Technit/ues 

The  flow  sections  of  the  geometries  of  Fig.  1  were  of  nonconduebng  material.  The  jet  of  cold  air  (hereafter 
referred  to  as  fuel)  emerged  from  a  2  mm  orifice  in  the  fuel  receiver  of  geometry  la  and  was  oscillated  by  a 
needle-valve  mounted  on  a  3(X)  W  vibrator  whose  di^Iacement  varied  direeby  with  input  voltage  from  a  power 
amplifier  and  inversely  with  the  square  of  the  frequency.  The  frequency  of  oscillabon  of  fuel  was  limited  to 
around  200  Hz  since  higher  frequencies  implied  smi^I  amplitudes  of  oscillabon  of  the  fuel  jet  In  geometries  lb 
and  Ic,  the  2  mm  fuel  jet  entered  radially  into  a  tube  of  diameter  7.4  mm  with  a  straight  upstream  length 
exceeding  100  diameters  to  ensure  fully  developed  flow  of  secondary  air  with  mean  temperature  up  to  45  degC 
higher  than  that  of  the  fuel.  The  lube  legih  downstream  of  the  fuel  jet  was  straight  in  geometry  lb  and  had  a 
bend  of  radius  25  mm  in  geometry  2c.  In  geometry  Id  the  secon^y  air-fuel  jet  entered  the  main  duct  of 
diameter  80  mm  coaxially  and  in  geometry  le  radially  at  the  wall.  Mcasuremenu  were  also  made  in  geometries 
Id  and  le  with  a  40  mm  flame-holding  ^sk  downstream  of  the  inlet  of  the  secondary  air-fuel  mixture.  The 
mean  flow  rates  of  the  fuel,  secondary  and  main  air  of  up  to  0.015, 0.18  and  2.7  kghnin,  respeebvely,  measured 
by  calibrated  float-type  flow  meters,  were  associated  with  mass  flow  rauos  of  jet  to  cross-flow  between  5  and  15 
and  Reynolds  numben  up  to  25  <X)0  in  the  secondary  air  lube  and  40  (XX)  in  the  main  duct,  corresponding  to 
turbulent  flow  condibons  comparable  with  those  in  the  acbve  control  measurements  of  refs  (3.4]. 

Hot-wire  anemometry  was  used  to  quanbfy  mean  and  rms  velocibes  in  geometries  la  and  lb.  Spectrum 
analysis  of  bie  amplified  hot-wire  signal  yielded  the  amplitude  of  the  fluctuabng  component  of  velocity  at  bie 
imposed  frequency.  The  rms  of  the  fluctuabng  component  at  the  imposed  frequency  was  less  than  2  %  of  the 
mean  velocity  compared  with  turbulence  intensity  of  order  4  %  at  the  exits  of  geomeoies  lb  and  Ic  due  to  the 
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small  ratio  of  fuel  to  tccoadary  Mir  Mild  damping,  and  mean  velocities  were  measured  using  Pitot  tubes  in 
geometries  lb  to  le.  Temperature  was  mcasui^  usasg  a  calibrated  cold-wire  resistance  thetmomeier.  The  rms 
temperuturc  at  the  imposed  frequency  was  obtained  by  spectrum  analysis  of  the  temperature  signal  in  preference 
to  using  probability  density  distributions  obtained  by  ensemble  averaging  because  of  the  shorter  sampling  period 
and  better  precision  in  the  presence  of  large  amplitudes  of  incoherent  fluctuations.  The  ampliuide  was 
determined  by  subtracting  the  spectral  power  at  die  imposed  frequency  &om  that  in  the  absence  of  imposed 
oscillations  and  was  subject  to  experimental  uncertainty  within  5%  in  flow  vrangeroents  lb  and  Ic  and  around 
10  %  in  I  d  and  le  associated  with  smaller  amplitudes  and  incoherent  oscillations.  The  am(ditude  of  oscillations 
across  a  flow  section  is  quamified  by  the  normalised  spatial  average  of  the  rms. 

5*wJpU  9dA/roj 

where  U  is  the  mean  axial  velocity,  p  is  the  density,  9  is  the  rms  of  the  oscillation  in  temperature,  t,  normalised 
by  the  temperature  difference,  AT,  between  the  secondary  air  iikI  fuel  and  m^  is  the  mass  flow  rate  of  fuel. 

RepeaubJiy  was  within  the  limits  of  experimental  uncertainty  and  further  details  of  the  measuiemem  techniques 
and  experimental  precision  are  available  in  ref.  (5) 

S.  Results 

Velocity  measurements  at  the  exit  of  the  orifice  of  geometry  la  conTirmed  that  the  amplitude  of  fluctuation 
of  the  jet  velocity  was  directly  proportional  to  the  vibrator  input  voltage  and  inversely  proportional  to  the  square 
of  the  frequency.  The  fluctuating  velocity  signal  at  120  Hz  was  sinusoidal  for  a  travel  of  the  needle  valve  ot  up 
to  0.48  mm,  estimatrxl  on  the  basis  of  the  characteristics  of  the  vibrator,  at  which  the  velocity  fluctuated 
between  zero  and  twice  the  mean  value.  The  results  of  Figs  2  to  6  are  for  geometries  lb  to  Id  and  a  mean  fuel 
flow  rau:  of  O.OIS  kg/min  oscillated  at  120  Hz  with  an  amplitude  of  0.015  kgAniit,  and  secondary  and  main  air 
flow  rates  of  0. 1 5  and  1 .35  kg/min,  respectively,  with  tecoodapr  air  temperature  30  degC  higher  than  that  of  the 

fuel  and  main  air  at  IS  Some  measureroents  were  carried  out.  especially  at  small  amplitudes  of  oscillation, 
with  a  temperature  difference  of  45  degC. 

Figure  2  shows  typical  distributions  of  mean  velocity  and  mean  and  rms  of  the  passive  scalar  at  the 
imposed  frequency  in  the  exit  plane  of  geometry  lb  along  two  orthogonal  diameters  for  two  values  of 
downstream  tube  length.  The  asymmetry  of  the  exit  velocity  and  oscillatory  scalar  profiles  in  the  vertical  plane 
for  the  shorter  tube  are  due  to  the  radial  jet  in  cross  flow  and  the  more  developed  (low  in  the  longer  duct  shows 
greater  symmetry.  The  amplitude  of  the  oscillations  decreases  sharply  in  the  vicinity  of  the  wall  due  to  the 
larger  radial  gradieni  in  pt^  caused  by  the  boundary  layer.  The  overall  attenuation  of  the  oscillation  is, 
however,  small  as  may  be  seen  from  Table  I  since  the  phase  difference  between  the  mean  flow  at  the  point  of 
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injection  and  the  exit  of  the  longer  tube  is  less  than  half  a  cycle.  The  amplitude  of  the  oacillaiion  at  220  Hz 
was  ahotu  5%  less.  The  diffetence  in  amplitude  is,  however,  small  and  of  the  same  order  of  magninide  as  the 
experimental  uncertainty.  The  exit  velocity  proGle  for  geometry  Ic  in  Fig.  3  shows  the  effea  of  the  bead 

associated  with  a  centripetal  acceleration  of  the  order  of  10^  m/s^  and  the  attenuation  of  the  oscillation  by 
around  40  %  was  mainly  due  to  the  the  Urge  gradients  in  velocity  and  phase  caused  by  the  reversed  flow  region 
along  the  inner  side  of  the  bead.  The  influence  of  the  mptiuide  of  osdllaiary  inpm  and  mass  flow  ratio  of 
secondary  air  to  fuel  on  attenuation  was  also  small  for  the  range  of  flow  eonditions  examined. 

The  use  of  flow  smoothing  screens  upstream  of  the  secondary-air  fuel  jet  resulted  in  a  nearly  uniform  mean 
velocity  distribution  in  the  main  air  flow  in  the  exit  plane  of  the  secondary  air-fuel  jet  in  geometry  Id.  Figure  4 
shows  radial  profiles  at  three  stations  and  axial  profiles  ttf  the  mean  velocity  and  mean  and  rms  scalar.  The  rates 
of  decay  of  the  mean  jet  velocity  and  scalar  are  in  accord  with  those  for  nirbulent  jets  and  the  imposed 
oscillations  in  the  jet  were  too  small  to  influence  the  spreading  of  the  jet.  The  slight  asymmetries  in  the 
profiles  in  the  horizontal  planearedueio  thaiinthe  jet  issuing  frm  a  tube  with  an  upstream  bend.  The  decay 
of  the  rms  scalar  is  faster  than  that  of  the  mean  piopenies  due  to  the  gradient  in  phase,  in  addition  to  the 
diffusion  effects  which  cause  a  decay  of  mean  properties,  and  increases  with  distance  from  the  axis  of  the  jet. 
Although  the  normalised  rms  of  the  flucniatioo  at  x  >  2.4  D  is  only  30  %  less  than  that  at  the  exit  plane  of  the 
jet,  (he  peak  amplitude  is  about  m  eighth  of  that  a  the  jet  exit  (see  Table  1).  The  rms  scalar  for  a  larger  main 
air  flow  rate  of  2.7  kg/min  was  around  S  %  less  than  that  with  the  mass  flow  rale  of  I.3S  kg/min  due  to 
increased  diluiioo  of  the  oscillating  flow.  The  difference  is,  however,  less  than  the  experimental  uncertainty  of 
around  10  %. 

Figure  S  shows  radial  profiles  for  the  mean  velocity  and  the  mean  and  rms  scalar  u  x  s  i.s  D  for  geometry 
le,  and  the  normalised  average  of  the  rms  scalar  is  less  than  a  third  of  that  at  the  jet  exit  (see  Table  1).  The 
oscillations  were  insignificant  at  x  >  2.4  D.  Radial  injection  leads  to  faster  diffusion  of  the  imposed  oscillations 
than  with  axial  injection  due  to  the  agiution  of  the  main  flow  by  the  radial  jet.  The  attenuation  was  also  an 
order  of  mngnilude  greater  than  in  geometry  Ib  due  to  the  mean  v^ity  in  the  80  mm  dua  being  a  tenth  of  that 
in  geometry  lb  so  that  radial  gradients  in  phase  develop  over  shorter  axial  distances 

The  location  of  a  disk  of  diameter  40  mm  (Q.S  D)  m  a  distance  of  O.I  D  Cram  exit  of  the  secondary  air-fuel 
jet  of  geometry  Id  resulted  in  severe  damping  of  the  osctUaiions  due  to  impingetnem  and  oscrUations  could  not 
be  detected  0.9  D  from  the  jet.  Although  Table  1  Aows  that  the  spatial  average  of  the  oscillatjons  in  the  plane 
of  the  disk  located  0.9  D  downstream  of  the  jet  was  only  S  %  less  than  that  without  the  disk,  the  maximum 
ampUiude  of  oscillation  in  the  pime  of  the  disk  wis  only  a  third  of  that  at  the  tame  axial  location  in  the  absence 
of  the  disk  owing  to  the  spread  of  the  jet  in  the  presence  of  the  bluff-body.  Fig.  6  shows  that  the  oscillations 
decay  rapidly  with  distance  downstream  of  the  disk  owing  to  the  efleci  of  the  turbulent  wake  behind  the  disk  on 
the  transit  of  the  oscillatory  scalar  and  oscillations  could  not  be  detected  beyond  1.2  D  from  the  disk.  The 
influence  of  the  disk  was,  however,  smaller  in  geometry  le  since  the  smpliuide  of  oscillations  were  largest  close 
to  the  dua  wall  opposite  the  point  of  entry  of  the  jet  and  away  from  the  disk. 

4  Discussioti 

The  mean  scalar  distribution  measured  20  tube  diameters  downsucam  of  the  oscillatory  jet  (geometry  Ib) 
indicates  good  mixing  and  the  normalised  spatial  avenge  rms  value  of  around  OJ  is  close  to  the  expected 
maximum  of  0.3S  corresponding  to  a  sinusoidal  oscillatory  input  The  attenuation  was  small  due  to  (he 
shortnas  of  the  tube.  The  damping  of  the  oscillations  near  the  wall  of  the  secondary  air-fuel  tube  and  in  the 
shear  layer  of  the  jet  in  parallel  flow  was  due  to  the  gradient  in  phase  caused  by  the  radial  variation  in  mean 
velocity  so  that  a  higher  Creqency  would  be  astociaied  with  a  larger  gradient  Although  differciKes  of  around  S  % 
in  attenuation  between  12(1  and  220  Hz  are  conaitieni  with  the  above,  measurements  with  a  wider  range  of 
frequency  values  and  longer  lube  lengths  are  necessary  to  be  more  conclusive.  Slow  fluid  flow  also  leads  to  an 
increase  in  gradient  in  phase  and  impingement  and  regions  of  rccirculaiory  flow  are  associated  with  Urge 
gradienu  in  phase  and  severe  atifiuiaiino  ^  oscilUdcns. 

Active  control  depends  on  (he  oacillaiion  of  beat  release  out  of  phase  with  the  insiahiliiy  and  a  large  local 
amplitude  of  oscilUtion  of  the  fuel  concentration  near  the  flame  from  could  be  more  advanuigeous  than  a  more 
even  distribution  of  the  oscilUtion  of  fuel,  which  will  also  be  subject  to  greater  incoherence  in  phase.  As  slow 
fluid  flow,  impingmem,  flow  reversal  and  agitation  of  die  flow  are  atsocUind  with  the  decay  of  the  oscilUiions, 
parallel  injection  of  the  oscillatory  air-fuel  jet  or  the  fuel  close  to  the  burning  zone  is  likely  to  conserve  the 
amplitude  and  coherence  of  the  oscilUiory  input  and  enhance  control. 
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Statiot 

t 

Tinax 

y 

lb,  L=65mm 

x,mm 

1 

x/D 

0.044 

0.307 

L*l45nini 

I 

0.037 

0.295 

lc 

1 

0.020 

0.176 

Id 

70 

0.9 

0.008 

0.145 

190 

2.4 

0.003 

0.12 

270 

3.4 

0.001 

- 

Id  with  disk 

70 

0.9 

0.003 

0.14 

le 

120 

1.5 

0.004 

0.1 

le  with  disk 

120 

1.5 

0.004 

0.1 

Table  1:  Spatial  average  of  rms  scalar 


Fig.  2;  Prorilas  of  mean  valocity  and  mean  and  rma 
scalar,  gaomatry  1b.  (a)  La65mm.  (b)  L>145mm; 
Kaimm,  O-buA  maan  valoci(y>55mia. 


Fig.  3;  ProWaa  of  maan  vatoeijr  and  maan  and  rma 
acalar.  gaomatry  lc.  x-1.  UaSSmA. 


Fig.  4;  Profilaa  of  maan  vabcity  and  maan  and  rmi 
scalar,  gaomatry  Id,  (a)  x-70nim,  (b)  »»190  mm. 
(c)  *-270mm.  (d)  axial  profila;  0-4.3m/s. 


Fig.  5:  Profilas  of  maan  vabcity  an^msan  and  rma 
scalar,  gaomatry  la.  x«120  mm.  U«4.3m/s. 


Fig  6:  Longitudinal  profilas  of  maan  valocity  and 
maan  and  rms  scalar,  gaomatry  1  d  with  40  mm 
disk  70  mm  downstraam  ol  jat,  y«30mm,  *«0. 
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